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ABSTRACT 

Aims. We characterize individual and ensemble properties of X-ray flares from stars in the Cygnus OB2 and ONC star-forming regions. 
Methods. We analyzed X-ray lightcurves of 1003 Cygnus OB2 sources observed with Chandra for 100 ksec and of 1616 ONC sources detected 
in the "Chandra Orion Ultra-deep Project" 850 ksec observation. We employed a binning-free maximum likelihood method to segment the 
light-curves into intervals of constants signal and identified flares on the basis of both the amplitude and the time-derivative of the source 
luminosity. We then derived and compared the flare frequency and energy distribution of Cygnus OB2 and ONC sources. The effect of the 
length of the observation on these results was investigated by repeating the statistical analysis on five lOOksec-long segments extracted from 
| the ONC data. 

Results. We detected 147 and 954 flares from the Cygnus OB2 and ONC sources, respectively. The flares in Cygnus OB2 have decay times 
ranging from <0.5 to about 10 hours. The flare energy distributions of all considered flare samples are described at high energies well by 
a power law with index a=-(2.1±0.1). At low energies, the distributions flatten, probably because of detection incompleteness. We derived 
average flare frequencies as a function of flare energy. The flare frequency is seen to depend on the source's intrinsic X-ray luminosity, but its 
determination is affected by the length of the observation. The slope of the high-energy tail of the energy distribution is, however, affected little. 
A comparison of Cygnus OB2 and ONC sources, accounting for observational biases, shows that the two populations, known to have similar 
X-ray emission levels, have very similar flare activity. 

Conclusions. Studies of flare activity are only comparable if performed consistently and taking the observation length into account. Flaring 
activity does not vary appreciably between the age of the ONC (~1 Myr) and that of Cygnus OB2 (~2Myr). The slope of the distribution of 
flare energies is consistent with the micro-flare explanation of the coronal heating. 

Key words, stars: activity, corona, low-mass individual: Cygnus OB2, ONC - X-rays: stars 
On-line material: machine readable tables. 



1. Introduction 

Pre-main sequence (PMS) stars have high levels of X-ray 
activity with non-flaring X-ray luminosities (L x ) up to 10 31 
ergs/s, about two magnitude orde r of above those observ ed in 



most main sequence (MS) stars (Preibisc h et al.L 12005b . The 



X-ray activity in the PMS phase is commonly attributed to a 
"scaled up" solar-like corona formed by active regions. For 
non-accreting PMS stars, i.e. weak T-Tauri stars (WTTSs), the 
fraction of energy emitted in the X-rays energies with respect to 
the bolometric luminosity (Lboi ) is close to the saturation level, 
~10" 3 , observed for rapidly rotating MS stars, supporting the 
idea of a common physical mechanism acting both in the MS 
and the PMS phases. A more complex scenario is observed for 
PMS stars that are still undergoing mass accretion via mag- 
netically funneled inflows (classical T-Tauri stars - CTTSs): 
soft X-rays are also produced here in accretion shocks (e.g. 
Argiroffi et all 120071) . but coronal activit y appears to be some 



what reduced with respect to WTTSs (IFlaccomio et al .. 2003 



Preibisch et all 120051) . Regardless of their accretion proper- 
ties, all PMS stars show the high-amplitude rapid variabil- 
ity associated with violent ma gnetic reconnection flares (e.g. 
Feigelson & Montmerlell 19991) . 

X-ray variability over a wide range of time scales and 
scenarios is common to all magnetically active stars (e.g. 
Favata & Micelal 120031: lGudell2004lstassun et alll2006h . On 



long time scales, this includes rotational modulation of active 
regi ons, their emergence and evolution, and ma gnetic cycles 
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(e.g. Ma rino et all 120031: IFlaccomio et all l2005h . Most of the 
observed variations, however, have short time-scales (~ hours) 
and can be attributed to the small-scale flares triggered by 
magnetic reconnection events. Various authors have proposed 
that the overall X-ray emission observed in magnetically ac- 
tive stars i s the result of a large number of ov erlapping small 
flares (e.g. Dra ke et al. , 2000; ICaramazza et all 120071) . in anal- 
ogy with the micro-flare hea ting mechanism proposed for the 
solar corona (lHudsonl[l99ll) . The energy distribution of these 
events is supposedly described well by a power law, such that 
the large majority of flares release only small amounts of en- 
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ergy, so that we only detect their integrated and t i me-ay eraged 
X-ray emission dAudard et al.ll2000tlGudel et al.ll2003l) . 

In the past, statistical studies of X-ray variability in low- 
mass stars were hindered by the poor temporal coverage of 
EINSTEIN and /?Q,SArob servations, which were often short 
and very fragmented (e.g. iFuhrmeister & Schmitt , l2003h . and 
also by the limited spectral coverage of these telescopes. They 
were essentially limited to soft energies, while flares are char- 
acterized by harder emission. However, since the launch of the 
XMM-Newton and Chandra satellites, these limitations have 
been eased. In fact, two major projects relevant to the study 
of PMS X-ray activity have been performed recently with the 
Chandra and XMM-Newton telescopes, targeting two galactic 
SFRs: 

i- The Chandra Orion Ultra-deep Project (COUP) consisting 
in an ~850ksec long (~10 days) Chandra observation of the 
Orion Nebula cluster (ONC), performed over a time span of 13 
days. With a total of 1616 detected X-ray sources, the ONC is 
so far the best-studied SFR in X-rays. The COUP observation 
has been used for a number of vari ability studi e s, amo ng which: 
flare statistics on 'young suns' dWolk et all 120051). p hysical 
modeling of inten se X-ray flares dFavata et all 12005 ). rota- 
tional modulation (Flaccomi o et alll2005l) . X-ray variability of 
hot stars Ste lzer et al.l d2005[). and flaring f rom very low mass 



(0.1-0.3 Mm) stars dCaramazza et all 120071) . 
ii- The XMM-Newton Extended Survey of Taurus Molecular 
Clouds (XEST). In contrast to the COUP, XEST consists of 
several different pointings wit h roughly uniform continuous 



exposures of ~ 30-40 ksec dGiidel et all 12007). An X-ray 



Stelzer et al 



variability study has recently been presented by 
(2007), focusing on the statistics of X-ray flaring sources in the 
context of the coronal heating processes. 

Cygnus OB2 is a massive ~2 Myr old star-forming region 
whose stellar population has recently b een studied in the X- 
rays by lAlbacete Colombo et all ( 2007 ). This work, based on 
a 100 ksec long Chandra observation, focused essentially on 
the detection of X-ray point sources, finding 1003 from the 
analysis of their X-ray spectra and on their optical and near- 
IR characterization. Here we present an X-ray variability study 
of the young low-mass stars of the Cygnus OB2 region, aimed 
at giving a statistical characterization of variability and, ulti- 
mately, understanding its physical origin. Through a compar- 
ison with other regions, this study will also allow us to un- 
derstand whether the X-ray variability properties depend on 
the different cluster environment and physical characteristics. 
Unfortunately, statistical results obtained for one region often 
cannot be directly compared with those obtained for another 
SFR, because of the different sensitivity limits and temporal 
coverages of the respective X-ray observations. However, us- 
ing the long COUP observation for the ONC, here re-analyzed 
consistently with the Cygnus OB2 data, we assess the effect of 
the observation length on our results and are able to present a 
robust comparison of flare variability in the two regions. 

In § 2 we describe the methods used to detect and charac- 
terize variability and, in particular, our operative definition of 
flares. In § 3 we present the observed flare properties: decay 
times, energy distribution and frequency. In § 4 we extend our 
study to the COUP data, considering both the entire 850 ksec 



observation and 5 distinct 100 ksec segments. Finally, we dis- 
cuss our results in § 5 and draw our conclusions . 

2. Variability analysis 

The analysis presented in this paper was performed starting 
from event lists of the Cygnus OB2 and ONC Chandra sources. 
These we re appropriately extracted in the 0.5-8. keV energy 



range by lAlbacete Colombo et al.l (120071) and iGetman et al 



(2005) for the Cygnus OB2 and ONC observations, respec- 
tively, in both ca ses using the ACIS -EXTRACT package 



dBroos et all 120021) . The analysis of the COUP data was per- 



formed in two ways: i- considering the entire 850 ksec COUP 
observation, ii- by selecting in the observations 5 different 100 
ksec segments of continuous observation starting at 0, 180, 400, 
650, and 850 ksec from the beginning of the observatiorQ. This 
approach proved useful for understanding the biases due to dif- 
ferences in the observation lengths. The 850 ksec COUP ob- 
servation, for example, per mits detection of ex tremely long du- 
ration (-400 ksec) flares dFavata et all l2005h . while in a 100 



ksec observation like the one available for Cygnus OB2 , we 
can only detect a fraction of the flares longer than 30-100 ksec 
and completely miss those with exponential decay times longer 
than ~100 ksec. This results in a bias against the observation 
of energetic flares. 

We initially searched the time series of each 
Cygnus OB2 star for variability using the one-sided 



Kolmogorov-Smirnov (KS) test ( Press et al. . 1992h . This 



test compares the distribution of photon arrival times with 
what is expected for a constant source and gives the con- 
fidence, Pks, with which the hypothesis that the source is 
constant can be rejected. Sources with Pks >99.9% have 
been considered as definitively variable, while those with 
99.0%<Pks^99.9% are considered as probably variable. 
Sources with Pks<99.0 are not considered significantly vari- 
able. In as Cygnus OB2 we found 135 sources out of a total of 
1003 sources with P KS > 9 9%, o f which 86 have P KS >99.9% 



(Alb acete Colombo et all 120070 . In the ONC the figures for 



the entire COUP observations (1616 detected sources) are 977 
and 886 dGetman et all l2005h . while for each of the five 100 



ksec segments we find, on average, 358 and 264 sources with 
Pks >99% and Pks >99.9%, respectively. Last two figures 
clearly show how the number of sources for which the KS-test 
detects variability in a given relatively short (e.g. 100 ksec) 
observation is a lower limit to the total number of variable 
sources in the region. This is essentially due to the fact that 
most of the observed variability is in the form of flares, i.e. 
events that are shorter than our observatio n and with duty 



cycles that are instead considerably longer dWolk et all 12005) 
than 100 ksec. Moreover, small flares and/or other low-level 
variability may remain undetected because the sensitivity of 
the KS tests critically depends on the source's photon statistics. 

We intend to specifically study flare variability here. In 
the next section we therefore briefly describe statistical tools 
that, in addition to indicating variability, provide an objec- 



1 Julian times (UT): 2452648.41, 2452650.58, 2452653.03, 
2452655.92, and 2452658.27 
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tive description of the time behavior of the X-ray emission. 
This, in conjunction with an operative definition based on our 
a priori idea of flare as an impulsive event, allows their effi- 
cient and, most importantly, unbiased detection in our observed 
lightcurves. 

2.1. Maximum likelihood block analysis 

We make use of a maximum likelihood algorithm that, un- 
der the assumption of Poisson noise, splits a light curve into 
periods of "constant" signal, referred to as blocks. In con- 
trast to more conventional approaches, this method works di- 
rectly on the sequence of photon arrival times and does not 
require binning. The maxim um likelihoo d block (MLB) algo- 
rithm is described by IWolk et al.l d2005l) . It has two free pa- 
rameters: the minimum number of counts per block (N m j n ) 
and the confidence level (CL) used in dete cting variab i lity. In 
the analysis of the 850 ksec COUP data, IWolk et alJ d2005b 
use the MLB algorithm with N m j n =20, because of the high 
photon statistic of solar-mass COUP sources. Unfortunately, 
our Cygnus OB2 observation is 8.5 times shorter than that of 
COUP, and the distance to the region is ~4 times greater than 
for the ONC. Most of the Cygnus OB2 sources have less than 
40 photons, compelling us to adopt a different N m j n for the 
analysis. A careful comparison showed that, for sources with 
more than 40 detected photons, variability is detected consis- 
tently using the MLB algorithm with both N min =l (MLB 1) 
and N m j n =20 (MLB 20). In sources with lower statistics, how- 
ever, detection of variability is hampered by MLB 20 because 
more than 40 photons are needed to define two different blocks. 
Thus, in practice, the MLB 1 algorithm is more sensitive to 
small flares and is just as effective as MLB 20 for more en- 
ergetic ones. 

Following IWolk et alJ (120051) . we classify blocks into three 
different groups according to their emission level: i— blocks 
compatible with the characteristic count rate (R c harX defined 
as the most frequent count rate exhibited by the source; ii— 
elevated blocks, during which the flux is above the R c har but 
usually not associated with impulsive events; Hi- very elevated 
blocks, marked by significantly elevated flux levels, often as- 
sociated with impulsive events. We also make use of a measure 
of the time variation in the count rate, i.e. the derivative, here 
defined as the ratio of the difference between the count rates 
of two successive blocks (AR=Ri-Pv2) and the minimum of the 
temporal lengths of the two blocks, mz'«(Ati,At2): 



dR 
dt 



l*i -Ri\ 

min(Ati, AI2) 



(1) 



We then define a flare as a group of non-characteristic blocks, 
beginning with a block for which the ratio between the deriva- 
tive and characteristic level is above a threshold, here set to 
10 



-4.2 , 



-1 . 



dR 1 

dt /? c har 



> Ah = 10 



-4.2 



(2) 



The duration of the flare, tfl r , is defined as the total dura- 
tion of consecutive elevated and very elevated blocks be- 
longing to the flare. Note that the adopted value of D t h is 




n. 



Time (Ksec) 




Time (Ksec) 

Fig. 1. Top: light-cu rve of a flaring Cygnus OB2 X-r ay source, 
#600 in the list of lAlbacete Colombo et ail d2007l) . The his- 
togram shows the light-curve binned in 1800 sec bins. The 
dotted-line shows the blocks as computed with the MLB al- 
gorithm. The first and last blocks are compatible with the char- 
acteristic level; the second is classified as "very elevated"; the 
third and fourth as "elevated" (see text). The blocks identi- 
fied as belonging to the flare are indicated by a thick solid 
segment at the top. The dashed line indicates an exponential 
fit of decay phase. Bottom: ratio between the derivative of the 
count rate and the characteristic level 4? • calculated at the 
boundaries of each block. The solid line indicates the thresh- 
old (10~ 4 2 ) used for the detection of flares. Note: The observed 
sequential second flare is an example of the sustained heating 
and therm odynamics cooling proc esses competing simultane- 
ously (e.g. iFavata & Micelal 120031) . This results in longer de- 
cay times. MLB compute these cases as single flares. 



slig htly lower than the one used by both Wo lk et al.l (|2005) 
and Caramazz a et al. (2007), i.e. 10 4 s The new, slightly 
lower threshold ensures the detection of 13 small but evident 
flares in the Cygnus OB2 data that would otherwise remain un- 
detected. Our block classification scheme is illustrated in Fig. 1 
for a Cygnus OB2 source with an evident flare. The top panel 
shows the light curve, both with fixed binning and in the MLB 
representation, while the lower one shows the time derivative 
at the interface between blocks. 

All the results presented throughout this paper are for 
CL=99.9%, N min =l, and D th =10- 42 . 
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3. Flare properties of Cygnus OB2 sources 

We applied the MLB algorithm and our operative definition 
of flare to the 1003 X-ray sources in the Cygnus OB2 region. 
We detected a total of 147 flares occurring in 143 X-ray 
sources. This means that about 14.7% of the sources appear 
to flare during the 100 ksec Chandra observation. We found 
that 74 of the flaring sources are also definitively variable ac- 
cording to the KS-test (P ks > 99.9%), while an other 25 have 
99%<P K s<99.9%. For the 48 remaining sources, P K s<99%. 
Most of these sources (32/48~67%) have less than 40 photons. 
Given that the confidence level used for the MLB algorithm 
was 99.9%, this seems to be more efficient than the KS-test 
in detecting variability, especially in the low-photon statistical 
regime. 

In order to estimate the energy of each flare, Efl r , we mul- 
tiply the number of flare photons, iVph-fk (=2 (fi ~ Rchar)&T), 
by a counts-to-energy conversion factor, CF. This is computed 
from the total energy released by the source during the obser- 
vation, i.e. the mean L x times the exposure time, divided by the 
number of detected source photons, A^ p h- sr c: 



ph-flr 



M>h- 



ph-src U0 



L T 

L x (t)dt = — — PjVph-fc = CF-AU-iir -0) 



Nph-src 



For Cygnus OB2 s o urces L x was computed by 



Albacete Colombo et al.l d2007l) using a single count rate 
for the flux conversion factor. For ONC, although X-ray 
luminosities from X-ray spectral fits for individual sources 



are available dGetman 



ray sp 
i et all 



2005), we derived new L x values, 



for consistency, using the same approach as was used for 
Cygnus OB2 . Thus Efl r was obtained in the two cases by mul- 
tiplying A^ph-flr by the relative counts-to-energy CF: 7.52xl0 32 
ergs/ph for the ONC and 7.92xl0 33 ergs/ph for Cygnus OB2 . 

We thus ignored both the source-to-source and the flare-to- 
flare variability in X-ray spectra. Given the low photon statistic 
of the sources and flares, in particular of those in Cygnus OB2 , 
it is indeed not possible to derive source-specific and flare- 
specific CFs. The dependence of CF on plasma temperature 
is, however, expected to be small. 

3.1. Flare duration and decay times 

Most observed flares appear to decay exponentially. We thus 
approximate their light curves as L x (t)=Loe~'/ T , where Lo is 
the peak luminosity and r the exponential decay time. The pa- 
rameter r is one of most important observationally accessible 
since it carries information on the physical conditions of the 
flarin g plasma and on the geometry of its confining magnetic 



field (ISerio et al.Lll991blReale et all 2004). 



Unfortunately, t is not easily measured for most of our 
flares because the photon statistic is insufficient for a reliable 
exponential fit. In the exponential assumption, the measured 
energy of a flare, Eg r , which we estimate from the number of 
photons detected from the beginning of the rise phase (assumed 
instantaneous) to its detected duration, Tg,-, can be written as 







L e-" T dt 



ir T(1 - 



-Tnr/r 



5 



■ 




Fig. 2. We illustrate flare parameters. Observed flare corre- 
spond to source N x 172. Values of the parameters can be fol- 
lowed in Table[3] Y-axis was converted from the count rate by 
using CF=7.92xl0 33 erg/ph. The shadowed area corresponds 
to the flare energy. Dashed line refers to the exponential decay 
fit of the flare. Note: Label "tau" in the equation refers to expo- 
nential decay time t. For flares characterized by three or more 
blocks, we give fitted values in Table[3] 

The additive term e~ T|, "' T in Eq.|4]could be understood as a cor- 
rection factor of t in a integration limited by Tg, . Note that the 
quantity ^ approaches r for Tg r » r. Hereafter, we indicate 

as te/l- In Fig[2]we illustrate computed flare parameters for 
a typical flare. 

In Fig. [3]we plot te/l vs. the time duration of flares, tflr, both 
estimated from the MLB analysis. We note in particular that Lo, 
entering into the calculation of Te/l. is approximated here with 
the maximum luminosity of the blocked flare light curve, L max . 
For an impulsive rise followed by an exponential decay, this 
last quantity is by definition smaller than the true peak flare lu- 
minosity (i.e. L max <Lo), but typically greater than -0.5 Lo. The 
vertical error bars in Fig.[3]reflect the resulting uncertainty on 
t E / L : 0.5E fll -/L max <TE/L<Efl r /L max . We also plot the theoretical 
loci, calculated from EqH] for values of t ranging from 1 to 24 
hours. A comparison of these theoretical curves with the data 
points indicates that the observed flares span a wide range of 
decay times. 

We now wish to test the usefulness of r™^ s (sEfl r /L max ) as 
a measure of the flare's decay time. For this purpose we have 
created a large number of simulated light curves assuming the 
same flare energy distribution found for Cygnus OB2 sources 
(see !33.2t and for a set of t ranging from 2 to 16 hours. 
We then applied the MLB algorithm and our flare defini- 
tion to characterize flare properties in the same way as for 
Cygnus OB2 sources. In Fig. H] we present histograms of t™^ s 
and tfl r resulting for simulations with r equal to 2, 4, and 6 
hours. While tflr appears broadly distributed, the distribution of 
t e/l S f° r simulated flare values peak at the corresponding t in- 
put value. This suggests that t™^ s may indeed be used as an es- 
timator of t, with a typical, associated uncertainty ct(t)~0.3t. 
Columns 12 and 13 of Table 1 give the values of t™™ s for all 
the 147 flares detected on Cygnus OB2 sources. The median 



■ , , is 2.9 hours and 90% of the values are in the [0.5,9] hour 
' range. For 27 flares for which the light curve comprised more 
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Table 1. X-ray properties of flaring sources in the Cygnus OB2 region. The complete version is available in the electronic edition 
of A&A 



N x 


Identification 


Albacete Colombo et al. (2007) 


KS-test 






MLB analysis 






# 


2MASS J+ 


Mass 


Phot. 


E x 


log(L x ) 


log(P K s) 


Rdmr 


T fc 


logdo 


Efli 


_meas 

'e/l 


rat 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


1 





— 


67 


3.23 


30.80 


-4.00 


0.158 


2.21 


32.10 


35.62 


0.91 


0.61 


6 


20322768+4113169 


1.59 


216 


2.61 


31.32 


-4.00 


0.958 


8.18 


31.98 


35.99 


2.87 


2.64 


8 


20322913+4114012 


0.48 


33 


2.62 


30.26 


-0.59 


0.300 


0.10 


31.46 


34.48 


0.30 


— 


20 


20323249+4113131 


2.27 


49 


2.64 


30.65 


-2.69 


0.145 


9.01 


31.01 


35.46 


7.70 


— 


28 


20323592+4112505 


0.94 


43 


2.04 


30.51 


-2.14 


0.245 


5.03 


31.15 


35.20 


3.06 


— 


33 


20323661+4122147 


6.83 


232 


1.93 


31.22 


-4.00 


0.985 


12.14 


31.74 


36.03 


5.40 


4.45 


37 


20323751+4117121 


0.18 


55 


2.68 


30.68 


-4.00 


0.104 


3.85 


31.44 


35.55 


3.59 


— 


39 


20323784+4122087 


1.41 


91 


2.20 


30.81 


-1.56 


0.463 


9.87 


31.15 


35.58 


7.43 


— 


43 


20323810+4112442 


1.19 


53 


2.74 


30.64 


-4.00 


0.246 


3.14 


31.38 


35.40 


2.91 


— 


51 


20323989+4113488 


0.58 


32 


1.97 


30.23 


-2.91 


0.127 


7.08 


30.90 


35.19 


5.46 


— 


52 


20324009+4110346 


7.97 


156 


2.69 


31.28 


-4.00 


0.533 


3.72 


32.04 


35.93 


2.14 


1.62 


54 


20324032+4122473 


1.16 


73 


2.68 


30.68 


-1.06 


0.570 


1.44 


31.50 


35.14 


1.21 




57 


20324063+4111113 


0.54 


30 


2.71 


31.05 


-2.31 


0.125 


8.68 


30.80 


35.19 


6.72 




74 


20324288+4119148 


1.52 


89 


1.86 


30.98 


-4.00 


0.238 


8.08 


31.44 


35.72 


5.25 


2.82 


81 


20324342+4111034 


0.60 


47 


2.20 


30.58 


-4.00 


0.180 


5.19 


31.15 


35.37 


4.55 




93 


20324492+4115203 


2.43 


45 


2.08 


30.99 


-1.52 


0.288 


6.22 


30.90 


35.13 


4.77 




125 


20324860+4108117 


1.27 


121 


2.23 


30.95 


-1.33 


1.165 


0.13 


31.68 


34.77 


0.34 




135 


20325011+4116270 


0.52 


13 


2.18 


26.58 


-1.48 


0.056 


3.76 


31.74 


34.80 


0.32 




145 






15 


2.49 


30.43 


-2.98 


0.046 


4.05 


30.85 


34.94 


3.43 




146 






27 


3.19 


30.32 


-0.59 


0.175 


2.21 


31.07 


34.92 


1.94 




159 


20325218+4118160 


2.83 


21 


1.79 


30.28 


-2.12 


0.112 


4.18 


30.83 


34.94 


3.59 




171 


20325351+4120195 


0.63 


11 


3.85 


31.99 


-0.89 


0.016 


4.72 


30.65 


34.89 


4.83 




172 


20325377+4115134 


2.39 


219 


2.42 


31.41 


-4.00 


0.924 


9.29 


31.80 


36.01 


4.47 


3.65 



Note: Rows refer to a source flare. Col. 1: source numbers according to Albacete Colombo et al (2007). Col. 2: name of 2MASS counterpart. 
Col. 3-6: source characteristics, i.e. mass of stars, "phot" the net detected photons, "E x " median photon energy, and "L x " X-ray luminosity 
[erg/s]. Col. 7: logarithm of P K s- Col. 8 count rate of the characteristic level in units 10~ 3 [s -1 ]. Col 9: flare duration. Col 10: logarithm of the 
unabsorbed L x at the flare peak was computed using CF. Col 11: logarithm of the flare energy (Eq.[3]l. Col 12: rL/E=Efl l /L pea ] [ . Col 13: fitted 
decay time of the flare. 



the median t™™ s and 
™l S f°i" tne whole 



than one block we also computed r from exponential fits to the 
block count rates. For these 27 sources 
Tfi t are both 2.6 hours. The distributions of t\ 
sample and for the 27 flares with exponential fits and that of 
Tfit for the latter sample are all statistically indistinguishable 
from each other using two-sample Kolmogorov-Smirnov tests 
(all giving null probabilities >16%). 

3.2. Flare energy distribution 



Stud ies of the sol ar corona (i.e. iLin et all [1984; 
Krucker & BenzL 1998 ) have found solid evidence of small- 
scale of continuous flaring. The distribution in energy of these 
flares was found to follow a power law: dN/dE oc E~ ff , where 
dN is the number of flares produced in a given time interval 
with a total energy (thermal and radiated) in the interval 
[E, E+dE], and a is the index of th e power-law distribution 
ilDatlowe et all 1 1974llHudsonl.il 99 ll) . 

The total energy released in flares is obtained by integra- 
tion: 



— EdE = k 



where E m i n and £ max are the minimum and maximum allowed 
flare energies, respectively, k is a normalization factor and in 
the last passage we assumed a>2 and £ max »£ m j n . For a>2, the 
integral diverges if E m \ n — >0, meaning that the energy released 
from micro flares (~10 27 -10 30 ergs) and nano flares (~10 - 
10 27 ergs) can contribute large amounts of energy to the to- 
tal emission an d eventually explain the whole coronal output 
dHudsonLfl99ll) . A non-zero value of E m [ n or a flattening of the 
flare energy distribution at energies is required to keep the to- 
tal emitted energy finite. The cumulative distribution of flare 
energies also follows a power law: 



E x - a dE 



~E 

a -2 " 



-(ot-2) 



■(5) 



N(> E ni 



r) <* | 

Je„, 



E~ a dE oc -E &1 



where j3 



(6) 



In the following, we compare the observed cumulative dis- 
tribution of flare energies with the above theoretical distribu- 
tion in order to constrain the slope, a. More precisely, we con- 
sider the observed frequency of flares per source vs. the min- 
imum flare energy, i.e. N(> £fl r ) divided by the number of 
sources in our sample and by the duration of the observatiorQ. 



2 We are using the ergodic principle here, assuming that all the 
sources in our sample share the same distribution of flare energies 
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Fig. 4. Distributions of the detected flare duration ta r (dashed lines) and t™£ s (solid lines) for simulated light curves constructed 
with flares having three different decay times r: 2,4, and 6 hours. Note that, while the distributions of t™^ s peak at the simulation 
input values, tft- appears broadly distributed towards higher values. 




T„ [hours] 

Fig. 3. TE/L(=Efl r /Lo) vs. Tf| r , the measured duration of flares. 
Filled circles indicate the observationally derived values, 
Eflr/Lmax, for flares detected on Cygnus OB2 stars; vertical 
error-bars the allowed range of t E / L considering the observa- 
tional biases on L : [0.5Eg r /L max ,Eg r /L max ]. Curves show the 
theoretical loci for exponentially decaying flares (Eq. HJ with 
different decay times, t, as labeled on the right-hand side (in 
hours). Note: Dashed line represents the limit case for which 
Tfl r » t (Eq. 2|, and thus the exponential correction factor 
(g-Tiir/r-) van i sn es (see text). 



Note however that, even if the simple power law model is cor- 
rect, the observed distributions will flatten at low energies due 
to an observational bias; i.e., small events are missed by our 



flare detection process. Above a certain threshold energy, E cut , 
however, our list of detected flares will probably be complete. 

Other than the sensitivity bias, a simple cumulative dis- 
tribution of flare energies can also be biased by the effect of 
overlapping flares. More specifically, the "effective" time avail- 
able for the identification of flares with a given energy is re- 
duced by the presence of larger ones, usually known as dead- 
time correction. This effect is generally small, given the low 
frequency of detectable flares. We correct for this bias, how- 
ever in the derivation of the frequency of flares by assum- 
ing that the time available for detect flares with a given en- 
ergy is red uced by the sum of the durations of more energetic 
flares (see, lAudard et all l2000t IStelzer et all I2007I) . Figure |5] 
(upper panel) shows, for the Cygnus OB2 X-ray sources, the 
dead-time-corrected cumulative distribution of flare energies, 
in units of flares per source per ksec. 

We use the max imum likelihood method described by 
Crawford et al. ( 1970h to determine E cut , i.e. the minimum en- 



(although they do not need to share the same rate of flaring). This is 
necessary in order to improve the statistics, given the low frequency 
of flares for a single source. 



ergy above which the distribution is compatible with a power 
law, and ft, the best-fit slope of the power l aw. The details of th e 
statistical analysis can also be followed in lStelzer et alj (12007b . 

In summary, we have computed, as a function of Efl r , 
the maximum-likelihood slope, yS(Efl r ), of the cumulative dis- 
tribution of flares with energies higher than Eg,. Alongside 
yS(Efl r ) we also used the Kolmogorov-Smirnov test to computed 
PKs(Efl r ), i.e. the probability that a power law with index y8(Efl r ) 
is indeed compatible with the observed distribution. The two 
functions are shown in the middle and bottom panels of Fig. [5] 
We chose E cut =10 35 1 ergs s _1 with respect to the maximum of 
Pks (-30%). This leads in an index 1.1 ±0.1. Note that fur- 
ther increasing E cut means /?(Efl r ) reaches a plateau, confirm- 
ing the compatibility of the observed flare energy distribution, 
above E cut , with a power law. 

In spite of our statistical analysis, our estimation of (5 and 
E cut from the observed energy flare distribution could still be bi- 
ased by the limited exposure time of the observation precluding 
the detection of a fraction of the longer flares. This bias could 
mainly affect the high-energy tail of the observed flare energy 
distribution of Cygnus OB2 sources. Unfortunately, with our 
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Fig. 5. Upper panel: Frequency of flares (flares per source per 
ksec) with energy above Efl r vs. Efl r , for Cygnus OB2 sources as 
derived from our 100 ksec Chandra observation. Middle panel: 
Maximum Likelihood estimate of the power-law index that bet- 
ter approximate the distribution for flare energies above Efl r . 
Three dashed lines, from top to bottom, ref fer to /lvalues 1.1, 
1.0, 0.9, respectively. Bottom panel: probability, according to 
the KS test, that a power law with best-fit index as given in the 
middle panel, is compatible with the observed distribution for 
energies above E cut . 



100 ksec observation, we are unable to quantify the degree of 
this incompleteness. 

In the extraordinarily long COUP observation (-850 ksec), 
however, this bias may be considered negligible. In $4] we 
present a re-analysis of the COUP data aimed at addressing this 
issue and at comparing the flare frequency of Cygnus OB2 stars 
with those of the ONC. First, however, we continue to discuss 
the issue of the frequency of flares. 



3.3. Flare frequency 

The low/intermediate-mass, X-ray-detected stellar pop- 
ulation within the Cygnus OB2 region covered by our 
Chandra observation counts , excluding 26 OB star s, 1003 - 
26=977 likely members dAlbacete Colombo et all 120071) . 



-2.0 - 



-2.5 



-3.5 - 



o 




-4.0 



35.0 35.5 
log (Efi r ) [erg] 



36.5 



Fig. 6. Flare frequency vs. minimum energy (as in Fig. [5]) for 
source sub samples with different minimum numbers of de- 
tected source photons (N p h), as indicated beside each curve. 
The dashed vertical line marks our flare detection detection 
completeness limit determined in § 13.21 Note: For sources with 
a similar N p h, differences in the energy of the flare are due to 
their respective exponential decay times (t). 



We exclude the OB stars from our statistics because their 
X-ray emission is mostly unrelated to magnetic activity. Our 
MLB analysis detected a total of 147 flares during the 100 
ksec observation. The average flare frequency (ffj r ) is thus 
approximately 1 flare in -664 ksec§ Considering only flares 
with energies above our completeness limit, E cut =10 35 1 erg, 
the flare frequency is reduced to ffl r =l/1320 ksec -1 . 

The flare frequency can, however, be expected to be a func- 
tion of source photon statistic, because of the combined effects 
of the limited sensitivity of our flare detection procedure and 
of the likely dependence of the flare intensity distribution on 
source brightness. Table [2] reports the flare frequencies calcu- 
lated for a sub sample of stars with different minimum num- 
bers of detected photons. Figure [6] also shows the full distri- 
butions of flare frequencies computed for different ranges of 
the source's photon statistics. Given the bias against the de- 
tection of low-energy flares, the flare frequencies for different 
sub samples are more directly comparable when considering 
only those flares with energies above our completeness limit, 
Eflr>E cut . This quantity is given in the last column of Table|2] 

Because of the strong dependence of the observed flare 
frequency on the sensitivity of the observation and on the 
intrinsic brightness of the X-ray sources, a direct compari- 



3 lAlbacete Colombo et all d2007ri indicate that up to -80 X-ray 
sources may be of extragalactic nature. Taking this contamination into 
account, the computed flare frequency should be slightly reduced to 
1/610 ksec- 1 . 
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Table 2. Flare frequency for stars with different photon statis- 
tics 



Min 


Nsrc 


Nflr 


f flr [ksec '] 


counts 


in sample 


in sample 


total 


comp.* 


1 


977 


147 


1/664 


1/1320 


10 


928 


143 


1/644 


1/1254 


20 


787 


140 


1/562 


1/1063 


30 


630 


132 


1/477 


1/863 


40 


531 


114 


1/465 


1/748 


50 


434 


102 


1/425 


1/629 


80 


265 


74 


1/358 


1/416 


120 


134 


44 


1/304 


1/352 


150 


78 


35 


1/222 


1/252 


180 


47 


33 


1/142 


1/162 


250 


13 


12 


1/108 


1/108 



The 26 X-ray sources identified with OB stars were excluded from the 
analysis of flare frequencies. Note: * indicates flare frequencies (ffl r J3 
with energies above our completeness limit (i.e. logEfl r >35.1 ergs). 



son with previous determinations for members of other SFRs 
is difficult. For example, the flare frequency (ff req ) computed 
for the PMS stars in the Taurus molecular clouds (TMC) by 
Stelzer et al.l (120071) is 1/200 ksec 1 , much higher than our fre- 
quency for the whole Cygnus OB2 sample. The difference may 
largely be explained as a sensitivity effect, given that the TMC 
is ~11 times closer than Cygnus OB2 and that the observa- 
tion were performed with a more sensitive instrument (XMM- 
Newton instead of Chandra ). Note that fo r a "complete" sam- 
ple of flares (i.e. those with Efl r >10 35 ergs~) IStelzer et al. (2007) 
find f flr ~ 1/770 kse c. During the COUP observation in the ONC 



Wolk et al. d2005l) observe 41 flares in a sample of 27 solar- 
mass stars implying ffl r ~ 1/650 ksec -1 . This figure reduces to 
fflr~l/l 150ksec~' if we only consider flares with energie s 
above the completeness limit estimated by lStelzer et alJ (2007). 
i.e. 10 35 3 ergs. 

Obviously, all these rates of flaring are not easily compared 
because of the different instruments, exposure times, sensitiv- 
ity, and the resulting completeness limits for the observations. 
The availability of the COUP data gives us the opportunity to 
re-analyze this observation consistently with Cygnus OB2 data, 
so as to minimize this differences. We give details of this anal- 
ysis in the next section. 

4. The analysis of the COUP data 

The COUP data set combines six nearly consecutive exposures 
of the ONC, spanning 13.2 days (1 140 ksec) with a total expo- 
sure time of 850 ksec (9.8 days). The observation was taken 
with the ACIS-I camera onboard Chandra . A complete de- 
scription of the COUP d ata analysis and sour ce detection pro- 
cedures can be found in iGetman et al.l (120051) . Here we make 
use of the 1616 source event files. 

To understand the bias introduced by the limited exposure 
time of our 100 ksec Cygnus OB2 observation, we present the 
analysis of the entire 850 ksec COUP observation here, as well 
as that of five different 100 ksec segments. As anticipated in 



$2] the 5x100 ksec segments (Seg. 1,2,3,4, and 5) were de- 
fined arbitrarily, but avoided time gaps in the COUP observa- 
tion. Figure [7] shows, as an example, the light curve of COUP 
source #1276. The horizontal segments labeled as "Seg. 1-5" 
indicate the time intervals for which we independently per- 
formed the flare analysis and we also indicate flares detected 
in the whole observation and in the 5 segments separately. 
Because the count rates obey Poison statistics, the maximum 
amplitude of fluctuations increases with exposure time. This 
implies that the statistical significance of a real signal (e.g. a 
flare) is higher when considering a shorter time interval. The 
MLB algorithm, applied with the same significance threshold 
(99.9% in our case) to the whole observation and to the shorter 
segments, can therefore yield different results. In particular, as 
exemplified in Fig. [7] faint flares may remain undetected in the 
former case (seg. 3), and consecutive flares may be detected as 
a single event (seg, 1 and 2). 

Inspection of the segmented light curves led to excluding 
three sources (COUP #9, #828 and #1462) from the following 
considerations, they lie close to the edges of the ACIS-I CCDs: 
because of the wobbling of Chandra , the light curves of these 
sources show high frequency variations at the wobbling period, 
leading to the spurious detection of a large number of flares. 

Figure [8] shows the rate of flaring vs. minimum flare en- 
ergy for the entire 850 ksec COUP observation, for each of the 
five independent 100 ksec COUP data segments, and for the 
combined flare population of the 5 segments. This curve, while 
maintaining the same properties with respect to observational 
biases, such the 5 curves for the single segments, improves the 
statistics and allows a more robust estimation of the mean flare 
frequency and of slope a. 

Results of the MLB analysis for the COUP data are pre- 
sented in Table[3] We report the following for each of the above 
data selections: N src , the total number of sources that show flare 
activity (row 1); Nfl r , the total number of detected flares (row 
2); l/ffl r , the inverse of the single-source flare frequency, i.e. 
the number of flares divided by the total observing time in ksec 
(row 3); E cut , the flare energy threshold above which flare de- 
tection is considered statistically complete (row 4); N^° mp , the 
number of flares with energies above E cut (row 5); l/ffj° mp , the 
inverse of the flare frequency for flares with energies above 
E C ut (row 6); a, the slope of the power law that characterizes 
the flare energy distribution above E cut (row 7); <x(a), the lcr 
uncertainty on a (row 8). Columns 2 to 6 give results for the 
analysis of the 5x100 ksec segments, column 7 the results for 
the summed 5x100 ksec segments, and column 8 the results for 
the analysis of the entire 850 ksec COUP observation. 

For the 1616 COUP sources observed for the whole 850 
ksec exposure, we found that detection of flares is complete for 
events with log(Efl r )>35.6, and the cumulative distribution fol- 
lows a power law with index q~2.05± 0.15, in agreement with 
results obtained bv lStelzer et al J (12007) for solar-mass stars (i.e. 
log(Efl r )>35.3 and a~ 1.9+0. 2). For the 100 ksec segments, we 
obtained a=\. 9+0.1, which is within ~lcr of the result ob- 
tained for the whole observation. Our analysis thus indicates 
that the slope of the flare energy distribution obtained from a 
100 ksec observation is statistically consistent with what is ob- 
tained from a much longer (850 ksec) observation. The flare 
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Fig. 7. Example of MLB flare detection for COUP source #1276. The result of the MLB algorithm applied to the whole obser- 
vation is overlaid on the binned light curve. The three thick horizontal segments on the top indicate flares detected using this 
representation. The four thick segments on the bottom instead indicate the flares detected by applying the MLB algorithm to the 
5 different 100 ksec time intervals, with latter indicated by the thinner segments labeled "Seg. 1-5". Note: time is given in hours 
from the beginning of the observation. 



Table 3. Statistics of X-ray flares in the ONC 



COUP 




5 x 100 ksec segments 




850 


Results 


Segl 


Seg2 


Seg3 


Seg4 


Seg5 


Total 


ksec 


N src 


125 


111 


110 


111 


127 


584 


640 




128 


112 


112 


115 


130 


601 


954 


1/fflr 


1253 


1432 


1382 


1395 


1233 


1334 


1429 


logE cut 


34.6 


34.6 


34.5 


34.6 


34.5 


34.7 


35.6 


•jcomp 


75 


63 


79 


62 


89 


313 


246 


i/C p 


2138 


2546 


2030 


2587 


1802 


2562 


5542 


a 


1.9 


1.9 


1.85 


1.9 


1.9 


1.9 


2.05 


cr(a) 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.15 



Notes: l/f flr is in ksec, E cut in erg. COUP sources #9, #828, and #1462 
were excluded from the analysis (see text). 



frequencies obtained from the shorter exposures are, however, 
significantly reduced. An inspection of the results of the flare 
detection processes reveals that many flares are missed in the 
short exposures either because the segment does not include the 
phases characterized by a high time derivative (the rise phase 
and the beginning of the decay) or, in a considerable number 
of cases, because the characteristic level is not observed and/or 
correctly determined. In some other cases, moreover, the flare 
is detected, but a large part of the flare falls outside the segment, 
thus underestimating its energy. 

4.1. Comparing the ONC with the Cygnus OB2 region 

Due to the small difference in age between the ONC (lMyr) 
and Cygnus OB2 (2 Myr), we expect the X-ray properties of 
stars in the two regions to be similar. Albacete Colombo et al. 
(2007) indeed found similar average X-ray emission levels. For 
the comparison of flare properties to be meaningful, we must 
consider results obtained from observations of the same tem- 
poral length. For the ONC we therefore consider the results 



-3.0 




34 35 36 37 

log(E fIr ) [erg] 

Fig. 8. Flare rate vs. minimum flare energy for the ONC X-ray 
sources observed by COUP, as determined from flares detected 
in: i- each of the five selected 100 ksec segments of the COUP 
data (dashed lines), ii- all the five segments (thin solid line), 
Hi- the whole 850 ksec COUP exposure (thick solid line). Note 
that, as discussed in the text, more low-energy flares are de- 
tected in the 100 ksec segments with respect to the whole ob- 
servation, but a significant fraction of the longer/more energetic 
ones are missed. 

from the the analysis of the five 100 ksec segments. Moreover, 
given the dependence of the X-ray luminosity on stellar mass 
and the relation between source luminosity and flare frequency 
( j33.3t . a meaningful comparison between the two regions 
should only consider sources in a similar range of mass. We 
therefore computed the flare energy distribution for the 92 
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Fig. 9. Rate of flaring vs. minimum flare energy for: i- the 
Cygnus OB2 sources with mass > 1 M observed during our 
100 ksec observation (thick solid line), ii- all the COUP sources 
observed during the five 100 ksec segments of the COUP ob- 
servation (dashed line), Hi- COUP sources with mass > 1 M Q 
in the five 100 ksec segments (thin solid line). 



ONC stars in the catalog of iGetman et al.1 d2005l) with masses 
> 1M , i.e. rough ly the lower mass limit o f X-ra y detected 



Cygnus OB2 stars dAlbacete Colombo et all 120071) . Figure 



compares this distribution with that of the Cygnus OB2 sources 
with mass >1M Q , both normalized to yield the average single- 
source rate of flaring. 

We analyzed the flare energy distribution for ONC sources 
with M > 1M using the same statistical procedure as described 
in §3. We find that, for Efl r >10 35 2 ergs, the distributions for 
Cygnus OB2 and ONC sources with masses higher than 1M Q 
(see thick and thin lines in Fig. [9] respectively) are compatible 
with a power law with index /?= 1 . 1 ±0. 1 (i.e. a~2. 1+0.1). We 
conclude that the flare statistics of the Cygnus OB2 and ONC 
stars are very similar. 

5. Summary and conclusions 

We conducted a systematic analysis of flare variabil- 
ity in the soft X-ray band of the young stars in the 
Cygnus OB2 star-forming region. For this purpose we analyzed 
the lightcurves of 1003 X-ray sources detected in a 100 ksec 
Chandra observation. For comparison we also used the same 
method to analysis of lightcurves of the 1616 X-ray sources 
detected in the ONC with the 850 ksec observation of the 
Chandra Orion Ultra-deep Project. To compare the results for 
the two regions, avoiding biases due to the different exposure 
times, we also analyzed, independently, five 100 ksec long seg- 
ments of the COUP observation. 

Flares were detected using the MLB algorithm. A to- 
tal of 147 flares were detected in the lightcurves of 143 



Cygnus OB2 sources. In the ONC, 954 flares were detected on 
640 sources in the analysis of the whole 850 ksec observation, 
while 601 flares were detected on 584 source in the analysis of 
the five 100 ksec segments. 

For each flare we estimated the emitted energy, Efl r , and 
the peak luminosity, L pea k. Backed by the analysis of exten- 
sive sets of simulated lightcurves, we suggest that the ratio 
Efli/Lpeak can be considered a reasonable estimate of the flare 
decay-time, t. This is particularly useful for weak flares, for 
which a fit to the light curve is not possible. The 147 flares in 
our Cygnus OB2 sample have a wide range of decay times from 
-0.5 hours to ~10 hours, with a median of 2.9 hours. For the 27 
flares that are bright enough, we fitted the decay phase with an 
exponential finding consistent results: a median t of 2.6 hours, 
and a distribution compatible with that of the whole sample. 

We find that Cygnus OB2 and ONC flare energy distri- 
butions display high-energy tails described by a power law 
(dN/dEccE~"). For energies below a given E cut , the distribu- 
tions flatten as a result of the incomplete detection of faint 
flares. For Cygnus OB2 sources, we obtain E cut =10 35 1 ergs and 
a-2. 1 +0. 1 . This slope agrees with the rang e of values found in 



previous studies of solar and stellar flares (IGiidel et all 12003 



Stel zer et all |2007|) and gives support to the micro-flare hy- 



pothesis for the explanation of the observe d X-ray emission 
and for the heating of corona dHudsoni 1 1 99 1 ). 

The average frequency of flares detected on any 
given Cygnus OB2 source during our 100 ksec observation is 
~l/664ksec~ 1 . It reduces to -1/1320 ksec -1 when considering 
only flares with Efl 1 >E cut =10 35 1 ergs, i.e. those for which de- 
tection is likely to be complete. It is, moreover, important to 
stress that our results, even for flares with energies higher than 
the "completeness limits", critically depend on the flare detec- 
tion method and on the duration of the X-ray observation. 

We investigated this point using the 850 ksec COUP data 
set, as well as for five distinct 100 ksec segments of the same 
observation. The frequencies as a function of minimum flare 
energy derived in the two cases are significantly different: the 
short exposures indeed hinder the detection of long and (usu- 
ally) energetic flares, often because it is not possible to cor- 
rectly determine the "characteristic level", which is instead 
clearly observed in the longer exposure. The "completeness 
limits", E cu t, above which the flare energies appear to follow 
power law distributions, are also different for the two cases: 
10 35 - 6 and 10 34 6 egrs for the 850 ksec and the 100 ksec ob- 
servations, respectively. The slopes of the power laws, a, are 
however notably similar, respectively 2.05+0.15 and 1.9+0.1, 
which are compatible with each other and with the slope of the 
Cygnus OB2 distribution within ~lcr. 

For the ONC stars we considered the analysis of the 
100 ksec segments. To compare stars with similar intrinsic 
X-ray luminosities, we restricted the sample to stars with 
masses higher than 1 M , i.e. roughly the completeness limit 
of our Cygnus OB2 observation. We find that the stars in the 
Cygnus OB2 and ONC regions have indistinguishable X-ray 
flare properties. 

Finally, we confirm that comparison of flare frequencies is 
only allowed if observational limitations and data analysis are 
performed in a single and homogeneous way. Contrary to this 
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result is just the determination of the slope of the power-law 
distribution, which is not critically influenced by the length of 
the observation. 
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